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Motivated by the recent scanning tunneling microscopy (STM) experiments [J. E. Hoffman et 
al., Science 297, 1148 (2002); K. McEIroy et al, Nature (to be published)], we investigate the real 
space local density of states (LDOS) induced by weak disorder in a d-wave superconductor. We 
first present the energy dependent LDOS images around a single weak defect at several energies, 
and then point out that the experimentally observed checkerboard pattern in the LDOS could be 
understood as a result of quasiparticle interferences by randomly distributed defects. It is also 
shown that the checkerboard pattern oriented along 45° to the Cu-O bonds at low energies would 
transform to that oriented parallel to the Cu-O bonds at higher energies. This result is consistent 
with the experiments. 

PACS number(s): 74.25.-q, 74.72.-h, 74.62.Dh 



Recently the energy-dependent modulations of the 
local density of states (LDOS) in high temperature 
superconductors has attracted a lot of experimental 
and theoretical attentions [1-6]. By employing the 
high-resolution Fourier-transform scanning tunneling mi- 
croscopy (STM), Hoffman et al. investigated the zero- 
field charge modulations in Bi2Sr2CaCu208+5 [1] de- 
duced from the checkerboard patterns in the real space 
LDOS. The charge modulation vectors near the origin of 
the momentum space were determined. They found that 
the period of the modulations depends on the energy and 
doping for energy below the maximum superconducting 
gap. With increasing energy (doping fixed) or doping 
(energy fixed), the LDOS modulation wave vectors ori- 
ented parallel to the Cu-O bonds become shorter while 
those along 45° to the Cu-O bonds become longer. They 
also observed that when energy increases, the charge 
modulation along 45" to the Cu-O bonds changes to that 
along the Cu-O bonds. Subsequently, McElroy et al. 
extended the previous measurements [1] to the second 
Brillouin zone and discovered the characteristic octet of 
quasiparticle states [2, 4]. 

A number of theoretical studies have devoted to the 
explanation of the STM experiments [3-6]. It has been 
proposed that the experimental phenomenon is due to 
the result of quasiparticle interference induced by disor- 
der [1, 2]. Following this idea, Wang and Lee calculated 
the Fourier component of the LDOS produced by an sin- 
gle impurity with a moderate strength on-site potential 
[3] . They obtained the LDOS images in momentum space 
and compared them with those in Rcf. [1], but didn't ex- 
amine the relations among the modulation wave vectors, 
dopings and the bias voltages. The present authors, on 
the other hand, investigated the effect of quasiparticle 
scatterings from a weak and extended impurity or de- 
fect [4]. The Fourier transform images of the LDOS and 
the relations among modulation wave vectors at different 
dopings and the bias voltages obtained by us [4] in the 



first Brillouin zone are consistent with the experimental 
observations [1, 2]. There exist also other works trying 
to understand the STM experiments [5, 6]. However, al- 
most all the previous studies [3,4,5] were restricted to the 
discussion of the Fourier transform of the LDOS due to 
a single defect in the first Brillouin zone. So far there 
have existed no studies of the real space LDOS images at 
different energies, and the origin of the checkerboard pat- 
terns observed in the experiments. Thus it is necessary 
to do these calculations, and to compare the obtained 
results directly with the LDOS images in the STM ex- 
periments [1, 2]. 

In this paper, we base on the approach of our previ- 
ous study [4] and examine the effect due to quasiparticle 
scattering from weak defects or impurities with both hop- 
ping and pairing modifications on the LDOS. The reason 
we choose weak defects as scatters is simply to elimi- 
nate the contribution from the resonant states. Using 
the T-matrix approach, We first calculate the energy de- 
pendent LDOS images due to a single extended defect 
and show how the pattern of the image changes as the 
energy varies. Then we demonstrate that the experimen- 
tally observed checkerboard pattern for the LDOS could 
be understood as a result of interferences among ran- 
domly distributed defects. Our results clearly indicate 
the dominant LDOS modulation along 45° to the Cu-O 
bonds at lower energy would transform to that oriented 
parallel to the Cu-O bonds at higher energy. This is also 
consistent with the STM experiments [1,2]. 

The Hamiltonian describing the scattering of quasipar- 
ticles from M impurities with local modifications of both 
hopping and pairing parameters in a d-wave supercon- 
ductor can be written as 

H = Hqcs + Hijnp, (1) 

where 

Hbcs = Xl^'^k - M)4<TCk<T + Au(4^cli^^ + c_kiCkT), 

kcr k 
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Here fj, is the chemical potential to be deter- 
mined by doping, = ti{coskx + cosky)/2 + 
t2coskxcosky + t3{cos2kx + cos2ky)/2 + ti(cos2kxCOsky + 
coskxCos2ky)/2 + tzcos2kxCos2ky, where ti_5 = 
-0.5951,0.1636,-0.0519,-0.1117,0.0510 (cV). The 
band parameters are taken from those of Norman et 
al. [7] for Bi2Sr2CaCu208+5, and the lattice constant 
a is set as a = 1. The order parameter away from the 
impurity is given by Ak = Ao(cosA;x — cosky)/2. 

Without loss of generality, at the impurity or defect 
site Ri, we assume an on-site potential consisting of a 
nonmagnetic part, Vsi, and a magnetic part, Vmi- The 
defect also induces a weak local modification in the hop- 
ping, Sti, to the nearest neighbor sites, and a suppres- 
sion of the superconductivity order parameter on the four 
bonds connected to the impurity site, SAn, and on the 
other twelve bonds connected to the nearest neighbor 
sites, dA2i. 

The Hamiltonian (1) with a single impurity has in 

fact been successfully applied by authors in Ref. 8 to 
explain the resonant STM spectra for Ni impurities in 
Bi2Sr2CaCu208-i-i. However, in the present case, no res- 
onances in LDOS have been observed in the recent STM 
experiments [1, 2]. So it is reasonable to assume that the 
on-site potentials (Vsi and Vmi) and the modifications 
in hopping and pairing parameters (Sti, 5 An and 5A2i) 
are all weak and have approximately the same order of 
magnitude. This model has been applied in Ref. 4 to 
explain the experimental observed Fourier transform of 
the LDOS at different energies and dopings in the first 
Brillouin zone [1, 2]. 

The Hamiltonian (1) can be solved by the standard Bo- 
goliubov transformation plus Green's function technique. 
When Sti,5Au,SA2i,Vsi and Vmi are all small, keeping 
the leading term in the T-matrix approach should be 
good approximation. In such an approximation, the res- 
onant states due to the impurity are eliminated and the 
LDOS change in real space due to these defects can be 
shown to have the following form 

M 

Sp{r,u:) = --^J2I1 E cos[(k-k')-(r-R,)] 

i=l k,k' v,v'=0,l 

x{[2^M(k,k') + 14i]a..' (k,k') 
+2[5AHB(k, k') + 5A2iC(k, k')]/3..' (k, k')} 

(2) 



where N is the number of sites in the lattice, A{\s., k') = 
coskx + cosky + cosk'x + cosk'y, i?(k, k') = coskx — cosky + 
cosk'x — cosk'y, (7(k, k') = cos{kx — 2k'x) — cos{kx 
■ — k'y) — cos{kx — k'^ + k'y) + cos{ky — k'^ — k'y) + 
cos{ky - k'y + k'^) - cos{ky - 2ky) + cos{2kx - k'^) + 
cos{kx + ky - k'y) - cos{kx + ky - k'^) - cos{2ky - k'y) + 
cos{kx - ky + k'y) - cos{kx - ky - k'x), Q;^^' (k, k') = 

Ck;/^k';/' ^ (-1)'^^'^ ^kiy^ki^+lCkV^kV + l, /^y;/' (k, k') = 

{-lY^^.u^^.u+l^l,, + (-i)'^'eLek',.'ek'.'+i, ci^iiuu^) = 

\/[iu}n — (— l)''i?k] is the bare Green's function, = 
V(ek - + A^, a. = [1 + (-l)''(ek - /x)/£;k]/2, and 
^0^1 = Ak/(2£:k). 

We note that Vmi is absent from Eq. (2) because there 
is no first order contribution from the magnetic poten- 
tial. Obviously, the total LDOS change Sp{r, uj) is a 
summation of those due to individual impurity. In the 
present study, we base our numerical calculation on a fi- 
nite lattice of 800 x 800 sites. For simplicity, we choose 
2Sti = Vsi = —25Au = — 4(5A2i, and assume that all 
these parameters are small such that the first order T- 
matrix approximation is valid. In our calculation, we also 
take the chemical potential ji = —0.1238 corresponding 
to the optical doping (15%) and introduce a finite life- 
time broadening 7 = 2 meV to the quasiparticle Green's 
function to smooth our data points by replacing uj + iO+ 
with UJ + i'y m Eq. (2). 

According to Eq. (2), We plot the images of the LDOS 
change due to a single defect (M = 1) located at the 
center of a 20 x 20 square lattice for different energies 
in Fig. 1. It is easy to see from the images at w = 
and -12meV that the LDOS modulations orient paral- 
lel to (±1,±1) directions (45° to Cu-0 bonds). When 
the energy becomes more negative at a; = —16,-20 
and — 25meV, the LDOS modulation clearly changes to 
(±1,0) or (0, ±1) (along the Cu-0 bonds) directions. At 
UJ =12 and 16meV, the LDOS modulations are strongly 
aligned along the directions of 45° to Cu-O bonds. For 
higher energies at w = 20 and 25meV, the modulations 
along the Cu-0 bonds begin to show up, and they may 
become dominant at higher energies. It is apparent that 
the LDOS images are asymmetric with respect to uj =0. 
We also note that with increasing energy \uj\, the region of 
the LDOS modulation along (±1, ±1) directions becomes 
smaller while that along (±1,0) or (0,±1) directions be- 
comes larger. 

In order to understand the LDOS modulations near a 
single impurity, we present the LDOS variations along 
(1, 0) and (1, 1) directions in Fig. 2. We can see that 
the LDOS on the impurity site has the maximum val- 
ues at energy uj < 12meV. At and above 12meV, the 
LDOS near the nearest neighboring site to the impurity 
have the maximum values. Far away the impurity (about 
20a), the LDOS modulations vanish. For w < — 12 meV, 



2 



the LDOS along (1,0) direction show strong and long dis- 
tance oscillations while those along (1,1) direction show 
rather week and short distance oscillations. This is the 
reason why the LDOS images shown in Fig.l in this en- 
ergy region have modulation vectors clearly along Cu-0 
bonds. At lo= -12, and 12meV, both modulations with 
approximately equal weight arc present in the LDOS. As 
a result, the LDOS images with modulations along the 
directions of 45° to Cu-0 bonds can also be seen in this 
region. For lj >12mcV, the strength of the modulation 
along (1,0) direction begins to outgrow that along (1,1) 
direction. We expect that dominant modulation should 
be along (1,0) direction, as the energy gets even higher. 
In fact there are total six nonequivalent charge modula- 
tion vectors [2,4], not all of them can be clearly identified 
in our real space studies. 

We have discussed the LDOS modulations due to a 
single impurity. In fact, there should be randomly dis- 
tributed defects in the experimental sample. In order to 
understand the effect of the disorder, we choose the de- 
fect concentration to be 1% and the defects are described 
by i?imp in Eq.(l). In our simulation, we first produce 
a random distribution for the defects and then calculate 
the LDOS changes in a 20 x 20 square lattice similar to 
that of the experiments [1, 2]. 

Graphs in Fig. 3 show the LDOS images at 9 different 
energies for such a distribution of defects. As a con- 
sequence of quasiparticle interference by these defects, 
checkerboard patterns in the LDOS show up in our nu- 
merical simulations. It is evident from Fig. 3 that the 
LDOS modulation or the orientation of the checkerboard 
pattern is parallel to the 45° direction from the Cu-0 
bonds at small |cj|, while at larger |a;|, the pattern tends 
to orient along the direction of the Cu-0 bonds. The 
changing of the orientation is particularly apparent when 
uj becomes more negative. This conclusion is consistent 
with the STM experiments [1, 2]. 

In siimmary, we have stiidicd the LDOS change in- 
duced by disordered and weak defects using Bogoliubov 
transformation and the the Green's function technique. 
The obtained LDOS images due to raudouily distributed 
defects exhibit checkerboard patterns which are simi- 
lar to those observed in the STM experiments [1, 2]. 
With increasing energy |a;|, the LDOS modulation along 
(±1, ±1) direction tends to orient itself in the (±1, 0) or 
(0, ±1) direction. This modulation transformation has 
also been seen in these experiments. Combining our pre- 
vious work [4], we conclude that the STM images [1, 2] 
can be qualitatively understood in the present theory. 
For the charge density wave order observed in Ref. [9] , it 
is believed to be the dimerization hopping and transverse 



pairing modulations [10, 11, 12]. However, the origin of 
such a static order needs to be further studied. 
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FIG. 1. The LDOS change 5p{r,u)) at different energy due 
to a single impurity at the center of a 20 x 20 square. 

FIG. 2. The LDOS change 5p{v,u)) versus the distance |r| 
to the single impurity along (1,0) and (1, 1) directions at dif- 
ferent energy. 

FIG. 3. The LDOS change 5p{v, w) at different energy due 
to the random impurities in a 20 x 20 square. 
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